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Abstract 

We perform event-driven molecular dynamics simulations of a system composed by uniaxial hard 
ellipsoids for different values of the aspect-ratio and packing fraction . We compare the molecular 
orientational-dependent structure factors previously calculated within the Percus-Yevick approx- 
imation with the numerical results. The agreement between theoretical and numerical results is 
rather satisfactory. We also show that, for specific orientational quantities, the molecular structure 
factors are sensitive to the particle shape and can be used to distinguish prolate from oblate ellip- 
soids. A first-order theoretical expansion around the spherical shape and a geometrical analysis of 
the configurations confirms and explains such an observation. 

PACS numbers: 61.20.Ja, 61.25.Em, 61.30.Cz, 61.20.Gy 
Keywords: 
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I. INTRODUCTION 



The structure of simple liquids can be characterized by the density- density correlator g(r) 
or by its Fourier transform, the static structure factor S(q). Strong efforts have been made 
to derive analytical tools for the calculation of these quantities, since g(r) or S(q) allow 
one to calculate several thermodynamical quantities, e.g. the equation of state. Among the 
most prominent theoretical approaches are the Percus-Yevick (PY), the hypernetted chain 
(HNC) and other more elaborate integral equations^. For models of simple liquids — when 
the number density p is not too large — integral equations provide predictions for g(r) and 
S(q) which well reproduce the "exact" results evaluated from experiments or simulations. 

For molecular liquids, structural information becomes more diverse due to the presence 
of the orientational degrees of freedom and of their interactions with the translational ones. 
Expansion of the angular dependent microscopic density with respect to spherical harmonics 
and Wigner functions for linear and arbitrary molecules, respectively, leads to a generaliza- 
tion of S(q) to tensorial correlators S\\>(q) (see section II. B). Several analytical approaches 
have been proposed to calculate these correlation functions. The simplest one performs 
"Ansatze" which relates S\\/(q), or the corresponding direct correlation function c\y(q), to 
S(q) or c(q), q = |q| of an appropriately chosen related simple liqui d 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 . Similar to 
this is the "Ansatz" for cw(q) based on the geometry of two molecule a 10 ' 11 ' 12 . The ther- 
modynamical perturbation theory is a systematic approach which uses a liquid system of 
isotropic particles as a reference and considers the deviation of the particles from sphericity 
as a perturbation. Different types of expansions exist, depending on the quantity which is 
expanded. For more details the reader may consult Ref.—. Finally, the integral equation 
theories have been extended to molecular liquidsi^. 

As compared to simple liquids, the quality of the integral equations for anisotropic par- 
ticles has been less intensively investigated. Furthermore, comparisons have been mainly 
restricted to hard ellipsoids of revolution and to the value of the coefficients gw m (r) and 
Cii' m {r) of the expansion into rotational invariants of, respectively, the pair distribution and 
direct correlation function in real space. The quality of PY and HNC-theory for a liquid 
of hard ellipsoids^ has been tested against molecular dynamic simulation data fifteen years 
ago^. Ref.— reports such a comparison together with earlier MC-resultsJ^, however, re- 
stricted to the center of mass correlator (?ooo( r )- Satisfactory agreement has been found for 
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both approximation schemes. PY-theory seems to be inferior to HNC-theory, because (i) 
the PY-correlators for oblate ellipsoids deviate stronger than the HNC-correlation from the 
corresponding MD-results and (ii) PY-theory does not yield an indication for an isotropic- 
nematic phase transition, in contrast to HNC-theory. The good quality of PY-theory for 
prolate ellipsoids has been confirme d 17 ! 18 . Particularly it has been demonstrated that rea- 
sonably good PY-predictions require that terms up to / max = 6 are taken into account. In 
variance with the observation of Ref.— it has been recently found that PY-theory also yields 
an isotropic-nematic phase transition 19 . 

Finally we mention that the direct correlation functions cw m (r) obtained from a MC- 
simulation 2 ^ have been compared with the "Ansatze" described in Refs 3 4 i 5 i 6 i 7 i 8 i 9 . Some of 
them exhibit a satisfactory agreement with the MC-results, particularly for large r-values. 
But in our opinion the agreement is less good than the one found from integral equation 
theor y 15 ' 17 ' 18 . 

The above short review on previous investigations of hard ellipsoids makes it obvious 
that the debate on the orientational structural properties is not quite settled. Therefore, we 
have applied a newly developed event-driven MD-algorithm 2 ^ in order to calculate the static 
molecular correlation functions of hard ellipsoids in q space. In contrast to previous work 
described above, we have chosen the tensorial correlators S\\>(q), A = (I, m) (see sect. II. B). 
They have the advantage that they can directly be deduced from neutron and synchrotron 
radiation scattering experiments, at least for I < 2. One of our main goals is to compare 
accurate numerical results with PY-predictions from Ref. [8], obtained for the gaussian 
overlap model 2 ^ and with a truncation at Z max = 8, a value larger than the one previously used 
in Ref.—. Besides this, we will interpret the peak structure of these correlators, particularly 
of the non-diagonal one, S2ooo(q)- We will demonstrate that the qualitative g-dependence 
of this nondiagonal correlator allows us to easily distinguish between oblate and prolate 
ellipsoids. The application of first order thermodynamic perturbation theory for 6*2000 (o) 
will provide support to this finding . Finally, we want to check how far the prediction of an 
isotropic-nematic phase transition 20 ' found from the growth of 62020 (q = 0) is reproduced by 
our MD-simulation. 

The outline of the manuscript is as follows. In the next section we will describe the model 
and the tensorial functions S\\> (q) . Section III contains results concerning: (i) a comparison 
for Sxy(q) fr° m PY-theory and MD-simulation, (ii) a geometrical interpretation of the q- 



4 



dependence of ^2000 (q) an d (iii) the first order perturbation theory for Sw(q). The final 
section IV contains our conclusions. 

II. METHODS 
A. Model 

We study a system composed of N = 256 or 2048 uniaxial hard ellipsoids, i.e. ellipsoids 
with revolution axis of length a and two other axes of identical length b in a volume V. 
The aspect ratio is defined as Xo = a/b, with < Xq < 00. The volume of each particle 
is 7rXo& 3 /6 and thus the packing fraction <fi = 7rXo& 3 p/6 (where p = N/V is the number 
density). We perform event driven molecular dynamics simulations, with periodic boundary 
conditions, at several values of 0.4 < <p < 0.51 and several 0.4 < X < 2.8 values. Distances 
are measured in units of the axes geometric mean I = \/ab 2 . Ellipsoids have mass m — 1 
and a spherically symmetric momentum of inertia, i.e. I x = I y = I z = 2mr 2 /5 with 
r = min{a, b}/2. 

The event-driven molecular dynamics simulation^ for a one-component hard-ellipsoids 
system (HES) used in this work is described in details in Ref: 23 . The prediction of events 
(collisions among molecules) is different from what has been proposed in the past. It relies 
on evaluations of distance between objects based on gradient descent and Newton-Raphson 
root finding algorithms. Such strategy in principle works for any objects whose surface is 
differentiable and hence it is not limited to ellipsoids. The efficiency of the algorithm is 
comparable to the algorithm recently proposed by Donev-Torquato-Stillinger— . 

To create the starting configuration, we generate a random configuration at very low 
<p and grow independently the particles up to the desired (i.e. we perform a quench 
in <p at fixed N,Xq). To equilibrate the systems, we simulate until the angular second 
Legendre polynomial correlation function, for the angle associated to the axis of symmetry, 
has decayed to zero. In addition we check that particles have diffused more than max{a, b}. 
Simulations last approximatively from 10 8 to 10 9 hard-ellipsoid collisions; between 1000 
and 5000 equally spaced configurations have been stored for the analysis. Fig. shows 
the studied state points together with the known thermodynamic lines. The majority of 
the studied points lies close to the equilibrium transition line to maximize the structural 
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signatures in the static correlation functions. 
B. Molecular structure factors 

A system of A" rigid molecules can be described by the positions of the centers Tj and 
the orientation (Euler angles) Qj of the j-th molecule. The microscopic density p (r, Q) = 
Y^j 5 (r — Yj) 8 (Q — Qj) can be expanded with respect to plane waves e 4qr and to Wigner 
matrices D l £ n (f2) 13 ,— . For molecules with a rotational symmetry axis, O = (#,</>), Wigner 
matrices reduce to spherical harmonics and the microscopic density can be expanded into 
tensorial modes 

where I takes integer values > and m runs between —I and I. The factor in front of the 
sum is for technical convenience. The molecular structure factors are defined as 

Siml'm' (q) = (P* m (q) Pl'm' (q)) 

that will in general depend both on the modulus and on the orientation of q. The symbol 
< .... > indicates an ensemble average. A convenient reference system is provided by the 
g-frame, where the direction of the z-axis is parallel to qi*- For molecules with a rotational 
symmetry axis, molecular structure factors in the g-frame become diagonal in so that 
SimVm' = 8mm' Simi'm- For ellipsoids, the Simi' m {q) have been calculated within the Percus- 
Yevick approximation^ and used as input of mode-coupling theory calculations to evaluate 
the glass transition lines^. 

III. RESULTS AND DISCUSSION 

A. Comparison between PY and simulation data for Si m Vm' 

Fig. El compares the numerical results and the PY predictions from Ref.— for 6*0000 , 5*2000 
and 5*2020 when X ps 1 (an almost hard-sphere case) for both oblate and prolate ellipsoids. 
In all cases, the PY predictions satisfactorily describe the numerical results. We note that 
5oooo resembles the typical shape of the HS fluid and is practically the same for oblate 
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(Xq = 0.9) and prolate X = 1.1) ellipsoids whereas S2020 is structureless for all q values (as 
expected since the studied points are far from the nematic phase). Nevertheless, information 
on the angular structure is contained in 5*2000- This function shows an interesting q behavior 
when comparing the prolate and oblate case. Indeed, it appears that the prolate S2000 (<?) 
has opposite sign as compared to the oblate one. We also notice that the location of the 
first peak of Soooo coincides with the location of one extremum also in the S , 2ooo(<?)- 

Fig- El shows data for the tensorial correlation functions for values deviating stronger from 
X = 1 namely X = 0.4 and X = 2.2 at = 0.4. Compared to the previous case, Soooo is 
less structured, while the opposite behavior is observed for both S2000 and S^o- Even in this 
case, prolate and oblate ellipsoids are easily distinguished from the q dependence of S 2 ooo : 
while for a prolate ellipsoids a maximum followed by a minimum is observed, the opposite 
behavior characterizes oblate ellipsoids. The location of the first peak of Soooo is now shifted 
as compared to the location of the S 20 oo(?) extremum. It is also interesting to observe that 
now a peak at q = is present in S2020, signaling the build-up of a finite nematic correlation 
length 27 on approaching the nematic transition. 

As a further case for comparing simulation results and theoretical predictions, Fig. 0] 
shows Si m ii m for (f) = 0.5 and X = 0.4 and 2.4. These state points are the closest points to 
the isotropic-nematic boundary for which PY predictions are available from ref.— . Even in 
this case, the PY results provide a satisfactory description of the q dependence. The growth 
of S 2 o2o (0) for q — > at high/low elongations (Fig 01 and Fig. 0J) signals the presence of 
a growing correlation length £2020 for nematic order. If £2020 would be of the order of the 
box-size, the results of our simulations would be affected by finite size effects. To test for the 
absence of such effects, we have simulated the two representative points <fi = 0.50, X = 0.40 
and 4> = 0.50, X = 2.80 for N = 2048 ellipsoids. As shown in Fig. El and Fig. El the data 
for the biggest size are consistent with the results for the smaller size at lower g's, providing 
evidence of the absence of finite size effects in the studied state point. Additionally, Fig. El 
and Fig. El depicts S2222(?) which has a pronounced peak at q — 0, and which is almost 
structureless for larger q. We have also calculated correlators with / and/or /' equal to 4. 
Those are not shown, because they are of less experimental relevance. Nevertheless, we 
mention that they satisfactory agree with PY- result. 
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B. Geometrical characterization of the q dependence of S2000 



An hand-waiving understanding of the anti-phase character of the oblate/prolate S2000 
function at same <fi can be obtained by a geometrical analysis of the configurations. For 
the case of prolate particles (X > 1), the location q* of the first minimum in 5*2000 (?) is 
located at q* ~ 2ii/b. At this small distance r* m b, the two particles must be almost 
parallel. Since q is along the ^-axis (g-frame) and the relative distance of the two particles 
has to be parallel to q in order to give a non-zero contribution to S2000, the polar angle 9 
of both ellipsoids is close to ir/2. This provides a negative contribution to S2000 because 
^2o(^?0) = (3 cos 2 9 — l)/2 « —1/2 (Fig. [7(a)) assumes its smallest possible value. Hence, 
around q*, where the majority of the pairs are parallel, S2000 will have a minimum. On the 
other hand, at the first maximum of S2000 (<? < <?*) correlation between pairs of ellipsoids 
with relative distance s* ~ max{a, b} > r* is sampled. As we are working in the g-frame, 
the pairs of ellipsoids at a distance s* that will contribute to the maximum at q* are those 
whose relative distance is parallel to q. If we take into account the excluded volume effects 
due to the ellipsoids at distance ~ r*, we see that the sterically favoured configurations are 
the ones shown in Fig. [7(b). For such pairs of ellipsoids, the total contribution to S2000 i n 
the g-frame is positive. Therefore if we find an extremum at a distance ~ s* for prolate 
ellipsoids we expect that one to be positive. 

The analysis of the configurations contributing to the peaks of S2000 f° r oblate ellipsoids 
is analogous. In the g-frame there is a shift of ±7r/2 in 9. Therefore, the sign of Y20 is 
inverted and, consequently the sign of S2000 is inverted with respect to the prolate case (see 



C. First order perturbation theory for Si m i> m i 

Specializing the results of Ref.~ to hard ellipsoids, it is possible to expand the radial 
distribution function g(r 12 , u 1; u 2 ) in a power series in e = X — 1: 



i=0 

where g n = 0(e n ), r\2 is the vector connecting two centers of the two ellipsoids and Uj are 
the unit vectors along the rotational symmetry axes. The first two terms of the expansion 



Fig. ED- 



00 




(1) 



s 



are 



0o(ri 2 ,ui,u 2 ) = 9Hs{ru) (2) 
#i(r 12 ,ui,u 2 ) = 

= -vQHsiruhi^, u 2 , e 12 )5(r 12 - a) 
-pa J d 3 r 3 g HS (r 12 ,r 13 ,r 23 ) ■ 

•[(7(ui,u 3 ;ei3)) U3 5(ri3 -a) + (3) 
+ (l(u 2 ,u 3 ,e 23 )} U[i 5(r 23 - a)} (4) 

were gHsij'12) and gHs( r '12, 1*13, 1*23) are the static 2-particle and 3-particle distribution func- 
tion of hard-spheres with effective diameter a depending on e, = r^e^-, p = N/V is the 
number density and 

(f(u)) u = — Jdvjd#Bin#f(u(#,(p)) . (5) 


The effective diameter a is defined as a = ((d(ui, u 2 , e 12 )) Ul ) U2 , where d(iii,u 2 ,ei 2 ) is the 
distance at contact of two ellipsoids with axis along u« and e i2 = V\ 2 jr\ 2 is the direction 
between the centers of the ellipsoids. The function 7 measures the "non-sphericity" of the 
potential and is defined self-consistently as d(iii,u 2 ,ei 2 ) = a [1 + 7(111, u 2 , ei 2 )] . At first 
order in e, 7 and a can be evaluated analytically from geometrical considerations 

a = b[l + e - + 0(, 2 )] (6) 

7(111, u 2 ,ei 2 ) = |[P 2 (ei2 -ui) +P 2 (ei 2 • u 2 )] +0(e 2 ) . (7) 

where Pi is the Legendre polynomial of order 

The general expression of the molecular structure factor Si m i>m in terms of <?(ri 2 , ui, u 2 ) 

is 

Sirrd'm(q) = 5 W + i 1 ' ~ l p Ait j d 3 r 12 ( ( g (n 2 , Ui , u 2 ) 
•e^ r -^(ui)^ m (u 2 ) ) U1 ) U2 (8) 



Using the properties of Legendre polynomials, we have that (7(11, w, e)) w = (e/3)P2(e , u), 
so we can rewrite g up to first order: 



g(r 12 , ui, u 2 ) = g H s(ri 2 ) - — {gHs{&)- [^(ui • e 12 ) + 
+P 2 (u 2 ■ eu)] 5(r 12 - a)+ 
+P J d 3 r 3 g HS (r 12 ,r 23 ,ri 3 )- 
[5(r 13 - a) P 2 (m ■ e 13 ) + 

6(r 23 - a) P 2 (u 2 ■ e 23 )] } (9) 

The zero-th order go is spherically symmetric and will contribute with a diagonal term 
Si mlm . The first order term contains functions of the form P 2 (u{ ■ e^) that can be recasted in 
terms of linear combinations of spherical harmonics Yj m (iij) with I = 0, 2; so the first order 
term will contribute only to Si m i' m with 1,1' = 0, 2 and I ^ V . 

In particular, in the g-frame the p-independent part of g\ contains only linear combina- 
tions of Yqo and F20; it is then possible to evaluate Si m i' m with the result 



Simi' m (q, e) = 5iirSf mlm (q) - e^-j=pa 3 g HS (a) ■ 
■ (x -1 — 2x~ 3 ^j sin a; + 2a; -2 cos 2 

" |yZm,2(A'm,00 + ^m,00^/'m,2o] 

+peFi m v m (<y) + 0(e 2 ) (10) 

where x = \q\a and J r imi' m (d') is a function of a that does not depend on e. Its calculation 
requires the knowledge of the static three-point correlator (7hs(^i2, ^13, ^23), which, however, 
is not known exactly. 

Therefore, 52000 shows a peak of opposite sign around q ~ a -1 as observed above when 
discussing the simulation data; in general, first order theory predicts 

D(q, e) = S 200 o(q, e) + S 200 o{q, -e) = (11) 

for e small enough. We check this last property by our simulations. Fig. El- (top) shows 
5 , 2ooo(q , ) e), S 2 ooo(q, — e) and D(q, e) for e = 0.05. The sum D(q, e) vanishes within the error, 
supporting the prediction of Eq. [TT] To gather a feeling of the range of validity of this 
prediction we show in Fig. El-(bottom) the average value (over all q < 50) of \D(q, e)| and 
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D(q,e) as a function of e 2 . The behavior of < D(q,e) > and < \D(q,e)\ > confirms the 
linear dependence on e 2 predicted from Eq. [TT]for e < 0.1. 

IV. CONCLUSIONS 

In this paper we have presented a comparison between PY predictions and MD simulation 
results for the tensorial correlators in a hard ellipsoids liquid, close to the phase coexistence 
lines. 

The major focus has been on the calculation of static molecular correlation functions. 
In contrast to earlier work we have not chosen the r-dependent coefficients gu> m (r) of an 
expansion with respect to rotational invariants but the tensorial correlators 5; m ;' m /(q) in 
g-space. Those have the advantage that they can be obtained from scattering experiments, 
at least for 1,1' < 2. The comparison of the correlators from MD-simulations with the 
corresponding ones from PY-theory^ is rather satisfactory for all correlators and all pairs of 
(Xo, (j)) we have studied. Accordingly, the good agreement between results from PY-theory 
and an earlier MD-simulation found in Ref.— is confirmed. An interesting observation made 
is the qualitative difference of the nondiagonal correlator 52ooo(<?) for oblate and prolate 
shape. Since the back transform to real space is a linear procedure, this qualitative different 
behavior should also exist for G^oooM which is related to the coefficient g 2 oo ( r )- Indeed, 
Fig. la for X = 2 and Fig. 4 for X = 1/3 from Ref.^ show that the first extremum of 
#200 ( r ) is a minimum and a maximum, respectively. 

The qualitative shape dependence of 5*2000(9) has been proven analytically. Using first 
order perturbation theory with respect to e = X — 1 we have shown that 

52000 (<?) I oblate ~ — 5200o(<?) I oblate 

We have not attempted to compare 52000(9) from this perturbational approach, with the 
corresponding result from our MD-simulation and PY-theory, because one needs the static 
three-point correlator for hard spheres as an input which is not known. 

It has recently been predicted 2 ^ that the time-dependent correlator S^oo^ t), which is a 
measure of the coupling between the center of mass (I = 0) and orientational ( "quadrupolar" 
part / = 2) motion, has an effect on the light scattering spectra. This effect has been found 
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experimentally^ and might offer the possibility to check how far spectra from light scattering 
experiments may allow to discriminate between oblate and prolate particles. 

Finally, we have checked the growth of nematic order. In contrast to the PY-result 
in Ref.— the authors of Ref.— have found that S2020 (<7 = 0) from PY-theory diverges at 
a critical volume fraction C (A O ) for X > 2 and X <• 0.5. At <p c an isotropic-nematic 
transition occurs. Our simulation reproduces the shape and the growth of the peak of 
5 , 202o(q') at low g's. We have also demonstrated that, for the chosen values of X and 0, no 
finite size effects influences this peak. We have not attempted to determine 4> C (X ) from our 
simulation since this type of analysis would request much larger system sizes than we are 
currently able to simulate. 

The PY tensorial correlators have been recently used as input in molecular mode-coupling 
theory to evaluate the glass lines in the (X — 0) plane. The theoretical calculations suggest 
the possibility of a new mechanism of slowing down of the dynamics driven by the increase 
of the nematic order. Therefore, the validity of the PY predictions, particularly for the peak 
of 5*2020 close to q = 0, presented in this work confirms that this new mechanism is not 
arising from a failure of the PY predictions, but it is a genuine prediction of the molecular 
mode-coupling approach. 
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FIG. 1: Hard ellipsoids phase diagram. The dashed line are the phase boundaries calculated bjj^. 
The open diamonds correspond to the points of the phase diagram which we compare with the PY 
results o^S. The X's correspond to the points of the phase diagram for which we have analized 
finite size effects 
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FIG. 2: Simi'm' (q) for Xq = 0.9 (top) and Xq = 1.1 (bottom) at <j> = 0.51, i.e. near the hard-sphere 
case Xq = 1. Symbols are simulation results, lines are PY predictions from Ref»^. 
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FIG. 3: Simi'm' (q) at (j) = 0.4 for Xq = 0.4 (top) and Xq = 2.2 (bottom). Symbols are simulation 
results, lines are PY predictions from Refj^. 
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FIG. 4: Simi'm' (q) at values of 4> near the nematic transition line, for values of X corresponding 
to oblate/prolate ellipsoids; Xq = 0.4, = 0.5 (top) and Xq = 2.4, <fi = 0.5 (bottom) 
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FIG. 5: Comparison of S2020 (q) (top) and S2222 (q) (bottom) for iV = 256 (full circles) and 
TV = 2048 (open circles) ellipsoids at <f> = 0.50, X = 0.40. 
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FIG. 6: Comparison of £2020 (q) (top) and £2222 (q) (bottom) for N = 256 (full circles) and 
N = 2048 (open circles) ellipsoids at (f> = 0.50, X = 2.8. 
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FIG. 7: (a) Sketch of configurations of two prolate ellipsoids with a center-to-center distance close 
to r*. The values of Y20 i n the q- frame is negative for the two prolate ellipsoids because their 
rotational symmetry axes are almost perpendicular to q. (b) Configurations of ellipsoids with a 
center-to-center distance corresponding to s*. The shaded region represents the volume excluded 
by particles at a distance ~ r*. The values of Y20 m the q- frame is positive because the rotational 
symmetry axes of both ellipsoids are almost parallel to q. The lines with the arrows placed on the 
ellipsoids centers represent the axes of symmetry of the ellipsoids. 
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FIG. 8: (a) Sketch of configurations of two oblate ellipsoids with a center-to-center distance 
close to r*. The values of Y20 i n the q- frame is positive for the two prolate ellipsoids because 
their rotational symmetry axes are almost parallel to q. (b) Configurations of ellipsoids with a 
center-to-center distance corresponding to s*. The shaded region represents the volume excluded 
by particles at a distance ~ r*. The values of Y20 in the q- frame is negative for prolate ellipsoids 
because the rotational symmetry axes of both ellipsoids are almost perpendicular to q. The lines 
with the arrows placed on the ellipsoids centers represent the axes of symmetry of the ellipsoids. 
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FIG. 9: Check of Eq^Jfor the case eft = 0.49. The upper panel shows £2000 (?) for Xq = 1 ± e and 
D(q,e) for e = 0.10. The lower panel shows the average of \D(q,e)\ (empty circles) over all q < 50 
from the simulation as a function of e 2 for several values of e. The behavior at low e is linear in e 2 
as predicted from first order perturbation theory. The presence of a constant term in < \D(q,e)\ > 
for e — > is due to the noise implicit in the measures. In fact, the values < D(q, e) > (empty 
squares) at small e are dominated by the noise and are scattered around zero. The e 2 regime seems 
to break down around e ~ 0.3 
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